Criterion for phase separation in one-dimensional driven systems 
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A general criterion for the existence of phase separation in driven density-conserving one- 
dimensional systems is proposed. It is suggested that phase separation is related to the size de- 
pendence of the steady-state currents of domains in the system. A quantitative criterion for the 
existence of phase separation is conjectured using a correspondence made between driven diffusive 
models and zero-range processes. The criterion is verified in all cases where analytical results are 
available, and predictions for other models are provided. 
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The existence of phase separation and spontaneous 
symmetry breaking in low-dimensional systems far from 
thermal equilibrium has been a subject of recent inter- 
est [0, ||. While it is well known that these phenomena 
do not take place in one dimension in thermal equilib- 
rium, several models of driven one dimensional systems 
with local dynamics have recently been demonstrated to 
exhibit both (|, ||, |j. Whether or not a given model 
exhibits phase separation is in many cases not a simple 
question to answer, and it may depend on numerical ev- 
idence which could be rather subtle. 

For example, in a recent 3-species model introduced 
by Arndt et al ||] (AHR), it has been suggested that one 
should expect two distinct phase separated states: one 
in which the three species are fully separated from each 
other (related to the phase separation observed by Evans 
et al |3) in a related model), and the other is a more sub- 
tle mixed state whose existence is supported by extensive 
numerical simulations of systems of finite length and by a 
mean- field treatment. Subsequently, an analytical analy- 
sis of the model has shown that the mixed state is in fact 
disordered, and that in order to see this in simulations 
one has to study extremely long systems (of the order of 
10 70 ), far beyond existing numerical capabilities [||. 

In another example introduced by Korniss et al a two 
lane extension of a 3-species driven system was stud- 
ied 0. It has been suggested that while for this model 
the one lane system does not exhibit phase separation [|| , 
this phenomenon does exist in the two lane model. The 
studies rely on numerical simulations of systems of length 
up to 10 4 . This result is rather surprising and not well 
understood. It may very well be the case that as for the 
AHR model, the two lane model does not actually ex- 
hibit phase separation in the thermodynamic limit, and 
that this could be seen only by studying extremely long 
systems. It would thus be of great importance to find 
other criteria, which could distinguish between models 
supporting phase separation from those which do not. 

In this Letter we introduce a simple general crite- 
rion for the existence of phase separation in density- 
conserving one-dimensional driven systems. Phase sep- 



aration is usually accompanied by a coarsening process 
in which small domains of, say, the high density phase 
coalesce, eventually leading to macroscopic phase sepa- 
ration. This process takes place as domains exchange 
particles through their currents. When smaller domains 
exchange particles with the environment with faster rates 
than larger domains, a coarsening process is expected, 
which may lead to phase separation. Our criterion quan- 
tifies this mechanism, and relates the existence of phase 
separation to the steady-state currents through which do- 
mains exchange particles. The criterion is readily appli- 
cable even in cases which cannot be decided by direct 
numerical simulations. 

In order to explicitly state the criterion we note that 
in many models which carry a non-zero current in the 
thermodynamic limit the current of a finite domain of 
size n takes the form J n = Joo(l + b/n) to leading order 
in 1/n. For b > the current of long domains is smaller 
than that of short ones, which leads to a tendency of the 
longer domains to grow at the expense of smaller ones. 
According to our criterion phase separation takes place 
at high-densities only for b > 2. Moreover if the current 
decays to its asymptotic value as J n = J 00 (l J rb/n CJ ) 1 the 
model is predicted to phase separate at any density for 
a < 1 while it is always homogeneous for a > 1. In some 
models Joo = 0, although, due to lack of detailed balance, 
the current J n of a finite system is non- vanishing. In this 
case the system is predicted to phase separate at any 
density. Models for which J n decays exponentially to 
zero with n have been analyzed in the past and indeed 
were shown to exhibit phase separation [|[ [|, ||. 

The results presented above emerge from a careful 
analysis of a zero-range process (ZRP) which could be 
viewed as a generic model for domain dynamics in one- 
dimension. To define this process we consider a one- 
dimensional lattice of M sites, or "boxes", with peri- 
odic boundary conditions. Particles, or "balls", are dis- 
tributed among the boxes with the box i occupied by rii 
balls. The dynamics is defined in the following way: a 
box i is chosen at random and a particle is removed from 
it and transferred to a left (right) neighbor with rates 



pw ni ((1 — p)w n .) where < p < 1. The rate w ni de- 
pends only on the number of balls in that box. The model 
may either be unbiased (p = 1/2), or biased (p ^ 1/2). 

In a grand canonical ensemble, namely an ensemble 
where the number of boxes M is fixed while the number 
of balls is allowed to fluctuate with their average number 
controlled by a fugacity z, the steady state weight of a 
configuration of the ZRP is known to be || 

M 

W ZRF ({n i }) = l[z n ^ nt . (1) 

1=1 

Here = Ylm=i ^/ w m for fc > 1 and JF = 1. In 
this ensemble boxes are statistically independent with 
a single-site occupation distribution function given by 
P{k) ~ z k T}~. Depending on the rates w n the model may 
or may not exhibit condensation in the thermodynamic 
limit, whereby the occupation number of one of the boxes 
becomes macroscopically large. Clearly the rate w n must 
be a decreasing function of n in order for larger blocks to 
be favored and to support condensation. It is known || 
that condensation occurs at any density when w n — > 
with n — > oo, or when it decreases to a non- vanishing 
asymptotic value Woo as (1 + b/n a ) with a < 1; no 
phase separation takes place for a > 1; for a = 1 phase 
separation takes place at high densities only for b > 2. 

This model may be used to gain physical insight into 
the dynamics of driven one-dimensional systems. Occu- 
pied boxes represent domains of the high density phase. 
The currents leaving domains are represented by the rates 
of the ZRP. This is done by identifying the rate w n asso- 
ciated with a box containing n balls with the currents J n 
leaving a domain of n particles. A bias in the currents 
to a certain direction may be incorporated through p as 
defined above. The existence of a box with a macroscopic 
occupation in the ZRP corresponds to phase separation 
in the driven model. 

In the following we consider several one-dimensional 
driven systems and study their domain dynamics by in- 
troducing a corresponding zero-range process. By ana- 
lyzing the ZRP, the existence of phase separation in the 
original model may be addressed. We begin by consider- 
ing the AHR model. We show that for this model the cor- 
responding ZRP yields its exact steady-state domain-size 
distribution. This ZRP does not exhibit phase separa- 
tion, in agreement with the results of Rajewsky et al [||. 
We then discuss the two-lane model and argue that it, 
too, does not exhibit phase separation contrary to results 
of numerical simulations of finite systems. 

The AHR model is a three-state model on a ring. Each 
site is either empty (0), occupied by a positive (+) or a 
negative (— ) particle. The model evolves by a random 
sequential dynamics in which a pair of nearest neighbor 
sites is chosen at random and exchanged with the rates: 
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FIG. 1: A typical configuration of the three- state model (bot- 
tom) and its corresponding configuration in the ZRP (top). 
Periodic boundary conditions are imposed on the two models. 

This dynamics conserves the number of particles of each 
type. As in most studies of this model we consider equal 
densities of positive and negative particles. Numerical 
simulations suggest that the model has three states a 
fully ordered state for q > 1, in which the system strongly 
phase separates into its three phases; a mixed state for 
Qc < Q < 1 m which the particles (both positive and neg- 
ative) in the system condense into a single high density 
phase separated from a low density gas-like phase; and a 
disordered state for q < q c where particles and vacancies 
are homogeneously distributed. On the other hand exact 
calculations within the grand canonical ensemble [|| show 
that the mixed state is in fact disordered, with a finite 
average length of the high density domains in the ther- 
modynamic limit. Thus the homogeneous and the mixed 
states constitute a single disordered state. The system 
therefore exhibits only two states, one fully phase sepa- 
rated for q > 1 and the other disordered for q < 1. 

It is known j|, |(J that for this model the steady-state 
weight, Wl(C), of a given microscopic configuration C is 

L 

W L (C) = Tr J] [z {S Ti (+)D + S^E) + S Ti(0) A] . (3) 

2=1 

Here L is the length of the lattice, tj — +, — , when site 
j is occupied by a +, — or 0, respectively and z is the 
fugacity which controls the average number of particles. 
The fugacity is the same for positive and negative parti- 
cles so their average number is equal. The matrices E 
and A satisfy 

DE - qED = D + E aDA = aAE = A (4) 

to give the correct steady- state weight. Explicit represen- 
tations which satisfy this algebra are known ||, ||. For 
our purpose it is sufficient to note that there exists a rep- 
resentation in which A = |0) (0| is a projection operator 
satisfying A 2 = A, where |0) is a state vector. 

We now consider the steady state of this model in an 
ensemble in which the number of vacancies M is held 
fixed while the number of particles, and thus also the 
length L of the lattice are allowed to fluctuate. We refer 
to this ensemble as a grand canonical. A typical configu- 
ration of the model consists of blocks of particles bounded 
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by vacancies. By a block we mean an uninterrupted se- 
quence of positive and negative particles, bounded be- 
tween two zeros. Let rii (i = 1, . . . , M) be the length of 
the ith block located to the left of the zth vacancy. The 
block lengths can take the values rii = 0, . . . , L — M and 
satisfy rii — L — M. The partial trace W({rii}) of all 
weights of microscopic configurations Wl(C) consistent 
with m, . . . , wm, takes the form 

M 

W({n i }) = J[z n >Z ni , (5) 

where Z^ = (0| (D + E) k |0) is the sum over all weights 
of microscopic configurations of a block of length k. Here 
we have used the above representation of the matrix A. 
Within the grand canonical ensemble the various do- 
mains are statistically independent with a domain size 
distribution P(k) ~ z k Z^. It is known that Z^, with 
the algebra (||), is the partition function (sum over all 
weights) of the partially asymmetric exclusion process 
(PASEP) on a one-dimensional lattice of k sites and open 
boundary conditions |ll|| . The boundary conditions 
are such that the positive particles are injected at rate a 
at the left end and are removed from the right end with 
the same rate. The dynamics in the bulk of the system is 
given by the same rates as in (Q). Moreover, the current 
in such a system of size k is given by J& = Zk-i/Zk so 
that Z k = Ylm=i l /Jm for k > 1. 

We now turn to define the corresponding ZRP. We con- 
sider a ZRP in which each box represents a vacancy in the 
AHR model. The box i is occupied by rii "balls" , which 
corresponds to the number of particles (positive and neg- 
ative) in the block to the left of the ith vacancy (see Fig. 
|l|). Since in the AHR model the positive and negative 
currents of a block of size n are equal to J n , we choose 
the hopping rates of the ZRP to be w n = 2 J n with a sym- 
metric hopping to the right and to the left. Here again 
J n is the current of the PASEP in an open system of size 
n. Thus for this ZRP one has = 2~ k Z]^. Hence the 
weight of a given block configuration ni, . . . , vim in the 
AHR and the corresponding configuration of the ZRP are 
the same (up to a 2 _fc prefactor which can be absorbed 
into the definition of the fugacity z) Q . The ZRP pro- 
vides a simple interpretation of the AHR steady-state 
dynamics in which blocks interact via exchange of parti- 
cles. Each block in the AHR behaves as a PASEP with 
open boundary conditions and equal injection and ejec- 
tion rates at steady state. Neighboring blocks exchange 
particles at a rate given by the PASEP current. 

We now use the ZRP to study the block size distribu- 
tion in the AHR model. To do this we use the asymptotic 
form of the current of the PASEP with q < 1, 

Jn = Joc(l + b/n + 0(l/n 2 )) , (6) 

with Joo = (1 - g)/4, b = 3/2 for 1 > q > 1 - 2a and 
b = -1 for q < 1 - 2a (H], Therefore, since b < 2 



in both regimes, no condensate can appear for q < 1, in 
agreement with Q. Moreover, using the form (|^) of J n 
it is easy to show that the block size distribution in the 
homogeneous phase is given by 

P(*)~^exp(-*/0 ; e = ^^. (7) 

It is evident that for b < 2 the average block size diverges 
as £ — > oo and the distribution function is valid for any 
density no matter how large. No phase transition takes 
place in this case. However, for b > 2 the average block 
size remains finite for £ — > oo necessitating the existence 
of a phase transition, which results in a macroscopic block 
at high densities. Note that at q = 1 — 2a the block-size 
distribution function changes in a non- analytic manner. 
This point was first noted for the case of q = in ||. 

This result yields interesting insight into the origin of 
the apparent transition seen in simulations whereby in 
a certain q interval the correlation length £ becomes ex- 
ceedingly large. It can be shown |1| that the q depen- 
dence of the correlation length £ is introduced by the 
higher order corrections (for example c/n 2 ) to the cur- 
rent @ . We have calculated the correlation length £ of a 
ZRP with rates w n = 1 + 3/2n + c/n 2 for a given density 
as a function of the parameter c (corresponding to chang- 
ing q) We find that £ exhibits a sharp increase of a 
few orders of magnitude over a narrow range of values of 
c. This reflects itself in large (but finite) blocks and an 
apparent phase separation in direct simulations. 

The physical picture emerging from this analysis offers 
a rather robust mechanism for phase separation, and we 
conjecture that it has a more general validity. We expect 
it to apply to other conserving driven models even though 
an exact correspondence to the ZRP may not be evident. 

We now demonstrate the use of the conjecture for the 
two-lane model introduced by Korniss et al j?J. This 
model is a generalization of (||) with q = to two lanes. 
Here in addition to the hopping process (||) within each 
lane, particles may hop to neighboring empty sites on the 
other lane with rate 7a, and to exchange with a neigh- 
boring particle on the other lane with rate 7. Numerical 
studies of the model have suggested that for large enough 
a the system phase separates . In this state the parti- 
cles condense into a single high density block. However, 
physical insight into the phenomenon is lacking. In par- 
ticular it is not understood why the two-lane model seems 
to exhibit phase separation while its singled-lane version 
does not ^. In order to apply the conjecture, the cur- 
rent J n of a block of size n is calculated for the two lane 
model. This is done by considering an open two lane 
system of length n with no vacancies where particles are 
injected and ejected at the boundaries with equal rates 
a. Setting 7 = 1 as in we show the results of nu- 
merical simulations for several values of a in Fig. |^. By 
comparing the results for the corrections to the current 
A n = (J n — Joo) /Joo with the line 2/n it is easy to see 
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FIG. 2: The finite-size corrections to the current A n = 
(J n — Joo)/Joo in the two- lane model with open boundary 
conditions, for different system sizes. Here 7 = 1. 



that in this case b < 2. In fact, the curve may be best 
fitted to b ~ 0.8. Our conjecture thus implies that phase 
separation does not take place in this model. Note that 
even if one tries to fit A n to l/n (T with d / 1 one finds 
that a > 1 yielding the same conclusion. We note that 
simulating the open systems and obtaining the asymp- 
totic behavior of the current involves a relatively modest 
numerical effort, as one only needs to simulate rather 
small systems. This should be compared with the huge 
systems which are needed in order to demonstrate the 
lack of phase separation in direct simulations. Indeed, 
the conjecture suggests that the apparent phase separa- 
tion found in numerical studies is due to simulations of 
systems much smaller than the typical domain size. 

We have studied another extension of model (||), 
whereby the hopping rates depend on both nearest and 
next-nearest neighbors. In these models we have shown 
that 6 = 3/2 with a = 1. The conjecture implies no 
phase separation, in contrast to our direct simulations. 

We now apply the conjecture to a class of models with 
vanishing J^. To this end, consider model ^ with q > 1. 
As before, we define a block as a sequence of particles, 
both positive and negative, bounded by two vacancies. 
The corresponding open system is the PASEP with q > 1, 
with open boundary conditions and equal injection and 
ejection rates, studied in JTT[ | . The current in such a 
block, as in the corresponding open system, is exponen- 
tially decreasing with the block size, J n ~ q~ n l 2 . This is 
easily understood, as the particles are moving against the 
bias when q > 1. Using the conjecture the system is ex- 
pected to phase separate for any density. Several similar 
models with exponentially decaying currents have been 
shown to exhibit a strongly phase separated state j|, ^, [| . 

We end by reiterating the assumptions involved in the 
conjecture, which relates phase separation to the currents 
J n of finite blocks. It is assumed that the current flowing 
through a block is given by its steady-state value and 
is independent of its neighboring blocks. This may be 
justified by the fact that the coarsening time of large 



domains is very long, and the domains have a chance to 
equilibrate long before they coarsen. 

Although the criterion introduced in this Letter has 
not been proved to hold in general, its underlying physi- 
cal mechanism for phase separation is rather robust, sug- 
gesting a broad applicability. It is of interest to test the 
conjecture on other models to check its general validity. 
In the models exhibiting phase separation we have stud- 
ied so far, the current vanishes. It would be inter- 
esting to find a one-dimensional driven model for which 
Joo > and b > 2. According to the conjecture such a 
model should exhibit a novel type of phase separation. 
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